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Abstract 

We study the propagation of neutrinos inside dense matter under the con- 
ditions prevaihng in a proto-neutron star. Equations of state obtained with 
different nuclear effective interactions (Skyrme type and Gogny type) are 
first discussed. It is found that for many interactions, spin and/or isospin 
instabiUties occur at densities larger than the saturation density of nu- 
clear matter. From this study we select two representative interactions, 
SLy230b and DIP. We calculate the response functions in pure neutron 
matter where nuclear correlations are described at the Hartree-Fock plus 
RPA level. These response functions allow us to evaluate neutrino mean 
free paths corresponding to neutral current processes. 



1. Introduction 



Neutrinos play a crucial role during the energy liberation phase that follows 
a super-nova collapse as they carry away most of the initial gravitational 
energy. In order to describe the evolution of the newly born neutron star 
it is important to have reliable estimates of the mean free path of neutri- 
nos in a medium with the typical characteristics of a proto-neutron star, 
namely with densities ranging between one and several times the value of 
the saturation density of symmetric nuclear matter, with temperatures of 
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around a few tens of MeV, and with a proton fraction which can represent 
up to 30% of the baryonic density. 

Wc report here on some calculations of the mean free path of neutrinos in 
pure neutron matter under various conditions of density and temperature. 
The scattering of neutrinos on neutrons is mediated by the neutral current 
of the electroweak interaction. In the non-relativistic limit and in the case 
of non-degenerate neutrinos, the mean free path of a neutrino with initial 
momentum ki is given by [1] 

1/A(ki,r) = -g^ydk3(^4(l + cose)5W(g,r) 

+4(3-cos0)5«(g,r)^ , (1) 

where T is the temperature, Gi? is the Fermi constant, cy (ca) the vec- 
tor (axial) coupling constant, ks the final neutrino momenta, q = ki — 
the transferred energy- momentum, and cos^ = ki • ks. The dynamical 
structure factors S^^\q,T) describe the response of neutron matter to ex- 
citations induced by neutrinos, and they contain the relevant information 
on the medium. The vector (axial) part of the neutral current gives rise 
to density (spin-density) fluctuations, corresponding to the S = (S = 1) 
spin channel. 

A microscopic framework based on effective nucleon-nuclcon interac- 
tions is used here to describe consistently both the Equation of State (EOS) 
and the dynamical structure factors of neutron matter. Then, the mean free 
path of neutrinos is calculated according to Eq. (1) for different values of 
the density and temperature of the nuclear medium. 

2. The EOS and the choice of effective interactions 

Nucleon-nucleon effective interactions are usually fitted so as to reproduce 
some selected properties of finite nuclei and symmetric nuclear matter at 
saturation. The study of proto-neutron stars requires the use of these inter- 
actions for nuclear densities well beyond the saturation density of symmet- 
ric nuclear matter, and we explain in this section the criteria for our choice 
among the plethora of effective interactions currently in use. In Figure 1 is 
represented the energy per particle in symmetric nuclear matter and pure 
neutron matter as a function of the density p, as resulting from a mean field 
calculation employing several Skyrme and Gogny interactions. For compar- 
ison, the variational results given in [2] are also plotted; such results are 
based on the two realistic two-body interactions AV14 and UV14, plus the 
three-body interaction UVII. The noticeable differences between these two 




Figure 1. Energy per particle of symmetric nuclear matter (left panel) and neutron mat- 
ter (right panel) as a function of density p for a selected set of Skyrme and Gogny effective 
interactions. The black squares and diamonds are the results of variational calculations [2] 
based on realistic two- and three-body interactions. 



upper bounds provide an estimate of the reliability of the calculations for 
densities beyond the saturation point. 

In the left panel of Fig. 1 one can see that the Skyrme forces SIII and 
SV overestimate E/A as compared with the variational results for values of 
the density larger than 0.25 fm~^. Such a discrepancy has to be related to 
the high value of the compression modulus provided by these interactions 
(~ 350 MeV). The remaining effective interactions used in this figure give 
more reasonable values for the incompressibility (around 240 MeV), and 
one can see that their predicted energies per particle are in reasonable 
agreement with the variational results in the range of densities displayed 
in the figure. In contrast, in neutron matter (right panel) only the effective 
interactions SLy230b and DIP lead to energies per particle in agreement 
with the variational results. As a matter of fact, both SLy230 and DIP have 
included the variational results of neutron matter into the adjustment of 
parameters. Thus, we select two representative forces, one of Skyrme type 
(SLy230b[3]) and one of Gogny type (D1P[4]). 

However, the agreement with the variational EOS is not the only cri- 
terium we have followed to select an effective interaction. We have also 
considered isospin and spin properties. The asymmetry energy is related to 
the second derivative of the energy per particle with respect to the asym- 
metry parameter / = [pp — Pn)/ P- For small values of /, one can write 

E{I)/A = E{{))/A + ei/2 + 0(/3) (2) 

where ei is the asymmetry energy coefficient. This coefficient must be pos- 
itive to guarantee the stability of symmetric nuclear matter, otherwise a 
small fluctuation of / would amplify leading to neutron matter as the more 
stable form of nuclear matter. In Figure2 is displayed the asymmetry en- 
ergy ei as a function of the density for the same interactions used in Figure 
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Figure 2. Asymmetry energy ei of nuclear matter as a function of density for a set of 
Skyrme and Gogny effective interactions. We also show the results of variational calcu- 
lations[2]. The black rectangle represents the empirical point including error bars. 




1. The value of ei as deduced from the semi-empirical mass formula (e ~ 32 
MeV at saturation, indicated in Figure 2 as a small black rectangle) is repro- 
duced by most of the effective interactions. Apart from the parametrization 
SLy230, the displayed Skyrme interactions lead to an isospin instability for 
densities a few times higher than the saturation density. The onset of such 
an instability is the value of the density at which ei becomes zero. 

We have also looked at the spin properties of the effective interactions. 
In particular, we have calculated the magnetic susceptibility of asymmetric 
nuclear matter [5] . In Figure 3 is plotted as a function of the proton fraction 
Xp = pp/ p the value of the ferromagnetic density, which is defined as the 
density at which a ferromagnetic instability appears. It can be seen that all 
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Figure 4- RPA dynamical structure factors S^^\q, T) for neutron matter as a function of 
transferred energy go- Calculations have been done at the saturation density of symmetric 
nuclear matter, for a transferred momentum q = fc_F/3 and for several temperatures. The 
interaction is SLy230b. 



Skyrme interactions predicts such an instability in the range of densities 
of relevance for the study of a proto- neutron star. Note that the case of 
pure neutron matter corresponds to Xp = 0, at the left of the figure. In 
contrast, Gogny-type interactions predict a ferromagnetic instability for 
densities which are well beyond the range of values where the mean field 
plus effective interactions can be reasonably used. 

In conclusion, we shall employ the interactions SLy230b and DIP for 
our calculations of the dynamical structure factors in neutron matter. 

3. Dynamical structure factors 

The detailed balance relation allows one to obtain the dynamical structure 
factors 5^"^^ from the imaginary part of the response function X^^Hq^T) in 
pure neutron matter as 

S^'Hq,T) = 1 — 9m x'^'\q,T) . (3) 

vr 1 - exp(-qo/T) 

We recall that q is a four- vector (qo, q), and this expression is valid for both 
positive and negative transferred energies qo- 

For Skyrme type interactions, the response function can be exactly cal- 
culated including the full RPA correlations (see [6] and references therein). 
The expressions for pure neutron matter are given in [6] (note that, in [6] 
Eq.(20) should be multiplied by a factor (1 - e-^^)/2). In Fi gure 4 are 
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displayed the dynamical structure factors for the density (S = 0) and spin- 
density (5=1) channels as a function of the transferred energy qo from the 
probe to the medium. At zero temperature, only positive values of the trans- 
ferred energy are accessible, which accounts for the heating of the medium 
by the probe. As the temperature increases, the structure factors grow in 
the negative transferred energy domain. This behaviour corresponds to the 
cooling of the medium by the probe. The pronounced zero sound in the 
spin-density channel at zero temperature is damped when the temperature 
increases. As a general rule, when the value of the density increases and 
approaches the ferromagnetic point mentioned in the previous section, a 
dramatic peak appears in the structure factor [6, 5]. 

We have also calculated the response function assuming several inter- 
esting approximations for the particle-hole effective interaction. For a free 
Fermi gas the response function is given by the well-known Lindhard func- 
tion. The effects of the interaction at the mean-field or Hartree-Fock level 
imply the use of an effective mass in the Lindhard function. This is exact 
for Skyrme-Hartree-Fock and represents a convenient approximation in the 
case of Gogny-Hartree-Fock[5]. The particle-hole interaction may be also 
approached by a Landau form, containing only I = multipole, or i = 0, 1 
multipoles. We stress that in all cases the particle-hole interaction has been 
consistently obtained from the particle-particle interaction used to describe 
the EOS. The full RPA response function has been exactly calculated for 
Skyrme interactions. For the Gogny interaction we have used the effective 
mass approximation for the Hartree-Fock response and the Landau approx- 
imation for the RPA response. 

4. Neutrino mean free paths 

The temperature dependence of the neutrino mean free path can be seen in 
Figure 5, where the calculated A has been represented in the case of a free 
Fermi gas, and in the mean field approximation using SLy230b and DIP 
interactions. We have chosen an incoming neutrino energy Ei = 60 MeV. 
Only at very low temperatures the behavior of A follows the expected law 
in [7] (indicated also in the figure). As the temperature increases the 
interaction effects, contained in the effective mass, become less important 
and one reaches the free Fermi gas value at very high temperatures. 

The effects of the residual interaction, i.e., of RPA correlations, can be 
seen in Figure 6, where we have represented the ratio of RPA to Hartree- 
Fock mean free paths as a function of the density. For SLy230b, we have 
calculated the response function using either the full particle-hole inter- 
action (labeled RPA in the figure) or its Landau approximation (labeled 
RPA-Landau in the figure). Only the latter approximation has been used 
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Figure 5. Mean free path in neutron matter as a function of the temperature. The curves 
labeled SLy230b and DIP have been obtained within the mean field approximation for 
the response function. The energy of the incident neutrino is Ei = 60 MeV. 
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Figure 6. The relative mean free path A^^^/A^^ as a function of the density for 
interactions SLy230b and DIP. 



for the interaction DIP. 

In the case of interaction SLy230b one can see that the effects of the 
residual interaction in Landau approximation (either / = or / = 0, 1) do 
not produce any significant difference as compared to the HF approxima- 
tion for densities around the saturation value. In the same region, the full 
residual interaction produces a sizeable increase of A. All three solutions 
predict a strong decrease of the mean free path when the density reaches 
the ferromagnetic density. This is due to the very strong axial response func- 
tion as the system comes closer to the ferromagnetic instability. In contrast, 
for interaction DIP a significant difference with respect to the HF results 
appears when the residual interaction is considered at the Landau level. 
Note that the differences in the results obtained using the Landau approx- 
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imations with / = or with / = 0, 1 are quite small. Finally, the force DIP 
also predicts a vanishing mean free path when the density increases, but 
the ferromagnetic instability happens at very high densities, well beyond 
the domain of applicability of the effective interaction approach. 

5. Conclusions 

In the domain of densities and protonic fractions that are relevant to the 
proto-neutron star evolution, the use of effective nucleon-nucleon interac- 
tions which arc originally designed for describing the neighbourhood of 
normal nuclear matter is a delicate issue. These interactions may fail to 
reproduce the EOS of neutron matter (SIII. SV, Dl, DIS), or they lead 
to isospin and/or spin instabilities at 2-3 times normal density (Skyrme 
forces) . 

Having this in mind, it is still possible to use some of the existing ef- 
fective interactions to calculate consistently the nuclear response functions 
and hence the neutrino mean free paths in neutron matter. The evolution of 
the response functions with increasing temperatures shows that neutrinos 
can efficiently mediate the cooling process of hot neutron matter. We find 
that RPA correlations have a significant effect on mean free paths. With 
Skyrme interactions, RPA correlations tend to increase mean free paths 
at densities around normal density whereas at larger densities they pre- 
dict very short mean free paths because of the onset of spin instabilities. 
Thus, these interactions do not favor a fast cooling process. On the other 
hand, the Gogny interaction predicts a quite different scenario where RPA 
correlations increase mean free paths in the whole range of densities (up 
to about 5po)- In this latter case where ferromagnetic instabilities are not 
present, the mean field and the RPA correlations make the nuclear matter 
more transparent compared to the free Fermi gas model. Thus, neutrinos 
would escape more quickly and the delayed mecanism will be less efficient. 
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